METHOD FOR TUNING PID CONTROLLERS 
APPLICABLE TO NONLINEAR SYSTEMS 



Field of the invention 

5 

The present invention relates to a method for designing 
PID controllers and, more particularly, to a method for 
tuning PID controllers which are applicable to nonlinear 
systems such as robot manipulators . 

10 

Background of the Invention 

A plurality of patents has disclosed gain tuning 
methods of PID controllers. These patents may be categorized 

15 into two types of methods; one is a gain tuning method of PID 
controllers using hardware equipments and the other using 
software algorithms . 

Among these patents, the second type of methods, i.e., 
gain tuning methods using software algorithms, to which the 

20 present invention pertains, can be exemplified as follows: U. 
S. Pat. No. 6,081,751 titled ''System and method for closed- 
loop autotuning of PID controllers", U. S. Pat. No. 5,971,579 
titled ''Unit and method for determining gains of a PID 
controller using a genetic algorithm", U. S. Pat. No. 

25 5,742,503 titled "Use of saturation relay feedback in PID 
controller tuning", U. S. Pat. No. 5,331,541 titled "PID 
control unit", U.S. Pat. No. 5,229,699 titled "Method and an 
apparatus for PID controller tuning", U. S. Pat. No. 
5,057,993 titled "Method and system for acquiring parameters 
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in process control", U. S. Pat. No. 4,754,391 titled ''Method 
of determining PID parameters and an autotuning controller 
using the method", and U. S. Pat. No. 4,466,054 titled 
''Improved proportional integral -derivative control apparatus". 
5 Based on online tuning, the above-mentioned patents are 

directed to methods for autotuning of PID controllers using 
algorithms which measure set point values and process 
variables to suggest suitable gain values. 

U. S. Pat. No. 6,081,751, as shown in Fig. lA, 

10 discloses a method for calculating new PID controller 
parameters either directly through the formulae associated 
with the Ziegler-Nichols reaction curve method or through the 
intermediate step of calculating an ultimate period and 
frequency from the time constant and dead time which are 

15 calculated from the period and amplitude of oscillation 
generated by a relay. 

U. S. Pat. No. 5,971,579 is directed to a method for 
determining gains of a PID controller utilizing a genetic 
algorithm unit shown in Fig. IB. 

20 U. S. Pat. No. 5,742,503 provides a method for 

autotuning parameters of a PID controller, wherein parameters 
of a transfer function are computed through two steps and 
precise parameters of the PID controller are detennined from 
the computed parameters of the transfer function. 

25 U. S. Pat. No. 5,331,541 discloses a PID control device 

which identifies the rise characteristics of a controlled 
system by a step response method on changing a reference, 
moves to PID control when idle time and slope successively 
obtained on the rise reach a predeteirmined value, and 
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computes PID control parameters based on the idle time and 
slope obtained up to that point . 

U. S. Pat. No. 5,22 9,699 suggests a method for tuning 
PID controllers, in which a proportional control gain is 
5 increased until a desired oscillation is obtained, an 
amplitude and period are measured from the oscillation, an 
ultimate gain and an ultimate period is calculated in 
accordance with the amplitude and period, and the parameters 
of the PID controller are set in dependent upon the ultimate 

10 gain and period. 

U. S. Pat. No. 5,057,993 introduces a method for 
acquiring parameters in the * process control . According to 
this method, a manipulated variable to which an 
identification signal from an identification signal generator 

15 is added is inputted to a process to produce a controlled 
variable output which is then sampled to obtain a dead time 
and a maximum gradient, and the initial values of PID control 
parameters are calculated on the basis of the dead time and 
the maximum gradient. The initial values of the PID control 

20 parameters and the like are set in the adaptation section. 
In the adaptation section, a pulse transfer function of the 
process is acquired, and PID control parameters are 
calculated from the acquired pulse transfer function by using 
a partial matching method in a frequency region. 

25 U. S. Pat. No. 4,754,391 discloses a method of 

determining the PID parameters for the PID controllers by 
monitoring a limit cycle generated in a controlled process to 
obtain characteristics of the process and determining optimum 
PID parameters to be used for succeeding process control on 
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the basis of the results of the limit cycle monitoring. 

U. S- Pat. No. 4,466,054 suggests a process control 
apparatus which comprises a nonlinear controller connected in 
parallel with a PID controller, a circuit for identifying a 
5 dynamic characteristic of a process and a circuit for 
determining gains of the PID controller according to the 
dynamic characteristic, which is illustrated in Fig. IC. 

All of the prior arts described above suggest methods 
for autotuning a PID controller to acquire a desired response 
10 during an online control process and realization of the 
methods in a hardware system. 

Such conventional methods, however, require additional 
systems, thereby making the structure of PID controllers to 
be complicated. Further, these methods suffer from the 
15 disadvantage of requiring a lot of effort for calculating 
necessary intermediate values . 

In addition, the object systems to which these 
conventional arts can be applied are limited to linear 
systems . 

20 

Summary of the Invention 

It is, therefore, an object of the present invention to 
provide a method for tuning a PID controller applicable to 
25 nonlinear MIMO systems represented in second order phase 
variable form without adding any supplementary devices to the 
PID controller. 

In accordance with an aspect of the present 
invention, there is provided a method for tuning a PID 
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controller, wherein the PID controller is comprised in a 
closed PID control loop system, the PID control loop 
receiving an input, the PID controller being coupled to an 
object system being controlled, wherein the object system 
5 outputs process variables which is supplied for comparison to 
the input, wherein a result of said comparison is supplied to 
the PID controller, the method comprising the steps of: 
inducing equivalent relationships between PID gains of the 
PID controller and parameters of time delay control (TDC) ; 

10 selecting a natural frecjuency vector and a damping ratio 
vector so as to acquire a desired error dynamics of the 
closed PID control loop system; selecting a sampling time of 
the closed PID control loop system; determining the 
parameters of TDC on the basis of the natural frequency 

15 vector, the damping ratio vector and a closed loop stability 
condition for TDC; and selecting PID gains of the PID 
controller on the basis of the equivalent relationships. 

Brief Description of the Drawings 

20 

The above and other objects, features and other 
advantages of the present invention will be more clearly 
understood from the following detailed description taken in 
conjunction with the accompanying drawings, in which: 
25 Fig. lA is a flow chart showing a conventional method 

for tuning a PID controller; 

Fig. IB is a block diagram showing another conventional 
method for tuning a PID controller; 

Fig. IC is a block diagram showing further another 
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conventional method for tuning a PID controller; 

Fig. 2 is a schematic view showing the overall 
configuration of a control system in accordance with an 
embodiment of the present invention; 
5 Fig. 3 is a block diagram showing a control object 

system and a PID controller in accordance with the present 
invention; 

Fig. 4A and 4B are flow charts showing the process of 
designing a PID controller; 
10 Fig. 5 is a perspective view showing a robot 

manipulator with six degrees of freedom which is an object 
system of a PID controller in accordance with the present 
invention; 

Fig. 6A is a graph showing response when a PID 
15 controller designed according to the present invention is 
applied to an object system subject to a step trajectory; 

Fig. 6B is a graph showing control input obtained when 
a PID controller designed according to the present invention 
is applied to an object system subject to a step trajectory; 
20 Fig. 7A is a graph showing response error when a PID 

controller designed according to the present invention is 
applied to an object system subject to a sinusoidal type 
trajectory; and 

Fig. 7B a graph showing control input when a PID 
25 controller designed according to the present invention is 
applied to an object system subject to a sinusoidal type 
trajectory. 
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Detailed Description of the preferred embodiments 



The application of the preferred embodiments of the 
present invention is best understood with reference to the 
5 accompanying drawings, wherein like reference numerals are 
used for like and corresponding parts, respectively. 

Fig. 2 shows a control system structure to which the 
present invention is applied. The control system includes a 
digital device (e.g., computer) 21, an I/O board 22, and an 
10 object system 23. A program for PID control runs in the 
digital device 21. 

The object system 23 is a nonlinear MIMO system in a 
second order phase variable form and may be represented by 
equation 1. 

15 

X -f A{xy x) = B{x, x)u ( 1 ) 

where, when the order of the system is n, x is an nXl vector 
of variables to be controlled, A{xyx) is an nXl nonlinear 

20 vector representing dynamic characteristics such as friction 
of the system, B{x,x) is an nXn matrix representing input 
distributions, and u is an nXl input vector of the object 
system. For instance, a robot manipulator may be represented 
by equation 1, in which case x is a vector showing the 

25 rotation angles for the individual axes of the robot 
manipulator. 

Referring to Fig. 3, the overall system of the present 
invention is shown in a block diagram where x^ is a desired 
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trajectory. 

First, the relationship between PID control and time 
delay control will be explained. Time delay control method 
is known as a robust control method for a sampled data system 
5 Based on the relationship, a gain tuning method in accordance 
with the present invention will be described in detail. 

The sampled data system includes both a discrete time 
system and a continuous time system. When the system 
represented by equation 1 is to be controlled, a digital 
10 device 21 such as a computer may be used as a controller. In 
this case, the digital device 21 is a discrete time system 
while the object system 23 is a continuous time system. 

As explained by K. Youcef-Toumi and S.-T. Wu in 
Input /output liberalization using time delay control -Trans . 
15 Of ASME, J. Dyn. Sys . , Meas . , Contr. , vol. 114, pp. 10-19 
(1992) , the time delay control (TDC) is a nonlinear robust 
control method. 

, When TDC is applied to the system represented by 
equation 1, it is expressed by the form of equation 2, below. 

20 

u{t) = u{t - ;L) + B'' (- x{t -A)-^x, (t) + K^e{t) + K^e{t)) ( 2 ) 

wherein e (t ) =xa (t) -x (t ) is an error vector, X is a time delay 
value, Ko and Kp are nXn constant diagonal matrix determining 
25 the overall closed loop system including the TDC and the 

object system to have desired error dynamics, and 5 ' is a 
parameter selected to satisfy equation 3, below 



I-BB''l < — ri r 1" (3) 

11.2 1 + [(1 + fi^rpjro + P2YPD K 



where the subscript i2 indicates an induced matrix 2 -norm and 
5 Yp=||K/.L ' yA^oL . Ypd=||Kp -k„^|.^ . 

The induced matrix norTti is defined as follows: As mXn matrix 
A of real elements defines a linear mapping y=Ax from into 
R", the induced p-norm of A is defined by 
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IK=sup\nF=maxlHI. (4) 



which for p=l,2,oo, is given by 



15 14 =maxEK| .ML =VnM^" . ML =max2;K| (5) 



where k^^^^J^ A) is the maximum eigenvalue of A^A. 

In equation 3, Pi and P2 are gains of L2 (see equation 9 
below) whose meanings are explained below. First, e(t) is 
20 defined by equation 6 

4/) = e(t) + Ki,e(t) + Kpe(t) (6) 

H and G is defined as H:d-^e and G:s\^e , respectively. 
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and defined as Ln^ norm of •(/) truncated at time T. 

Considering the operator H±:Si\->ei for each element of the 

vector £(t), the transfer function between e± and ei can be 
expressed by equation 7 : 

4i = *'W=^t/t^ (7) 

Likewise, considering the operator Gi:£i\—> , the 
transfer function between Si and can be expressed by 

equation 8: 

where koi and kpi are i*^*^ diagonal entries of Ko and Kp, 
respectively. Then, the gain of the transfer function L2 is 
defined by equation 9 : 

\\H,\l=max\h.{jo>} 
||G,||3=max|g,0fi>)| 

Further, ||//||, =||C7„||,., and = ||G„ ||„ . Also, {M„), =max\h,i/co} and 

(Mq)^i =max\gi(j£oi . Thus, pi and p2 are given by equation 10: 
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(10) 



For a sampled data system, the time t is expressed by 
the combination of the sampling time (At) of the control 
5 system with the step number (a=l , 2 , k) , that is, t=a-At. 
Accordingly, where TDC is used in the sampled data system, 
the time delay value X is the sampling time At of the control 
system. Further, x{t-X) is expressed by using the central 
differential method which is the error-lowest numerical 
10 differential methods while Xj{t) and ^(r)are expressed by using 
the backward differential methods, as shown in equation 11: 
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f\ xAk)-2x^(k-l) + Xj(k-2) , , 



At' 



At 



15 After siabstituting the parameters of equation 2 by- 
equation 11, rearrangement into the form of the PID 
controller makes equation 12 : 

V A^ At j 



In result, equation 12 is a TDC form of a sampled data 
system. 
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Generally, a PID controller 31 may be expressed in a 
form of equation 13 : 



(13) 



where K is an nXn constant diagonal proportional gain matrix, 
Td is an nXn constant diagonal matrix representing derivative 
times, and Ti is a constant diagonal matrix representing a 
reset or integral time. 

However, in order to show the relationship with TDC, a 
PID controller which includes a DC component vector will be 
considered, as illustrated in equation 14: 



where the DC component is a constant value determined by the 
initial error (e (0) , e (0) ) . 

In consequence, designing the PID controller 31 means 
selecting K, To, Tj and DC. 

Application of the Laplace transform to equation 14 
gives equation 15: 



Multiplication of s to both sides and rearrangement change 
equation 15 into equation 16: 




(14) 



U{s) = 1 + v + ~r/0^(5) + 



(15) 
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sU{s) = K(s + T^s^ + )e{s)-¥DC- sKT^eiO) (16) 

The inverse Laplace transform of equation 16 gives equation 
5 17: 



u{t) = K[e{t) + r^e(/) + r/'e(0)+ {DC + K{e{0) + T^eCO)) - u{Q))d{t) (17 



where e(0), e (0) and u(0) are values of e(t), e(t) and u(t) 
10 at t=to, respectively, and 8(t) is a Dirac delta function. 
In equation 17, DC is obtained by equation 18: 

DC = -/s:(e(0) + roe(0)) + tt(0) (18) 

15 where u (0) =Utdc (0) . 

Accordingly, equation (17) is simplified to equation 

19: 

«(0 = K[e{t) + T^m + Tr'eit)) (19) 

20 

Then, equation 19 may be rewritten to equation 20: 

m = K(ix, (0 - m) + {X, (0 - m) + t,' {x, (o - x(j))) (20) 

25 x(t) , 'x{t) , Xjit) , Xj(t) are obtained by the backward method, 

as shown in equation 21. 
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i(t) ^(A , x{k)-x{k-\) xik)-2x{k-l)+x{k-2) 

At At At^ 

(21) 

At At At^ 



In a sampled data system, a PID controller can be 
expressed as equation 22 : 

5 

u{k) = u{k-i)+Atu{k) (22) 

When equations 20 and 21 are substituted into equation 22, 
equation 23 is obtained: 

10 

u(k) = u{k - l) 
(23) 

Therefore, for a sampled data system, the PID 
15 controller can be expressed as equation 23. 

Comparison of equation 12 with equation 23 shows that 
the PID controller and the TDC take the same form in a 
sampled data system. Hence, Relationship 1 can be given 
between the gain of the PID controller and the parameter of 
20 the TDC: 



K 

Relationship 1: K= ^ T^=K^~' T.^K^ K'' 

At 
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Also, the DC component can be obtained from equation 

24 : 

5 DC = -K{e(0) + 7-^6(0))+ /i: ■ Af(e(0) + T,''e(0)) (24 ) 

In Relationship 1, At is the sampling time of the 
control system, Kd and Kp are parameters which define the 
error dynamics of the overall closed loop system and are 
10 determined by the intention of the designer. Thus, if only 

B ' is determined, all of the gains of the PID controller can 
be automatically determined , 

As mentioned above, 5 'is selected to satisfy equation 3 
and, in general, determined by the following two methods: 

15 (1) B^= a (a is a constant, I is a unit matrix) 

(2) 5 ' is a constant diagonal matrix (( B^ ) ii=cxi 
(i = l,...,n) ) . 

A detailed design process of the PID controller based 
on Relationship 1 will be explained with reference to Figs. 
20 4A and 4B. 

First, an MIMO type object system represented as a 
second order variable form is selected (Step SIOOO) . 
Relationships between parameters of PID and TDC are derived 
in connection with the object system (Step S2000) . 
25 Then, a desired error dynamics of the object system is 

determined (Step S3000) . The error dynamics is determined by 
the designer by selecting the natural frequency vector (co) 

15 



and the damping ratio vector (q) . and Kp are determined on 

the basis of the natural frequency vector and damping ratio 
vector selected in Step S3000 (S4000) . When i^^ terms of the 
natural frequency vector ©n and the damping ration vector q 
5 are (On± and qi, respectively, i*^^ diagonal entries of the 
matrixes Kb and Kp can be obtained through the relations 
kDi=2(;iCDni and kpi=C0ni^ . Afterwards, the sampling time of the 
control system is determined (Step S5000) . The sampling time 
is preferably set to be a small value, but influenced by the 
10 CPU speed of the digital device 210. 

In Steps S3100 and S3200, B ' is determined. By one of 

the two methods mentioned above, B * is set to satisfy 
equation 3. Finally, gains of the PID controller are 
selected on the basis of Relationship 1 and equation 24 in 
15 Step S3300. 

EXAMPLE 

An example will be described, in which the method of 
20 the present invention was applied to a virtual object system 
on computer. 

A Puma type robot manipulator having six degrees of 
freedom, like that shown in Fig. 5, was selected as the 
object system. This robot is a model depicted in the article 
25 "The explicit dynamic model and inert ial parameters of the 
PUMA 560 arm", IEEE Int. Conference on Robotics and 
Automations, pp. 510-518 (B. Armstrong, O. Khatib, and J. 
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Burdick (1986) ) . The dynamics of the robot manipulator is 
expressed by equation 25: 

M{e)e + v[9, e)+ g{g)-\- f(0,0)=t (25) 

5 

where M(0) is a 6x6 inertia matrix, 0 is a 6x1 vector 
representing the rotation angles with respect of six axes, 
v{0y0) is a vector representing Coriolis force and centrifugal 

force, G(9) is a gravity vector, f{0,0) is a 6x1 vector 

10 representing forces which are not included in system modeling 
such as frictional force or disturbance and t is a 6x1 torque 
vector applied on joints. Comparing equation 1 with equation 
25, one can get B (x, x) ->M (9) . 

Detailed design process will be described according to 

15 the process of Figs. 4A and 4B. 

For convenience, assume that poles of error dynamics of 
all six axes are located as multiple roots at 5. In this 
case, the natural frequency a)ni=5 (i=l,...,6), damping ratio 
<;i=l (i=l, 6) . Here, kDi=10 and kpi=25 (i=l,..-,6) in all of the 

20 six axes. The sampling time of the controller is At = 0.001 
sec. Puma 560 has an arm part and a wrist part. Since the 
two parts greatly differ in inertia due to the difference in 
mass and structure between themselves, two values (ai and aa) 
were used in which ai is for axes 1, 2 and 3 and a2 for axes 4, 

25 5 and 6. 
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||/-M-'Af|| < — f7 ^ T ).||/-Af-'Af|| <0.9732 



M = 



«1 -^xS 



0 

«2 -AxS 



->a, =1.8 flj =0.2 , 



|/-M-'m|| ^ = 0.7663 < 0.9732 . 
From Relationship 1 and equation 24, 



K, =K, =K, =^ = 18000,r„, =T^, =r^3 =-L- = 0.1,7), =7>, =7)3 =^ = 0.4 

=i:3 =^6 =^^ = 20oo,r^, = =r^, =-J- = o.i,r,, =r,, =r,, =^ = 0.4 



r 1 ^ 

it) = K, (e, (0 + I (o") Jcr + J^.e, (0 



Eighteen gains of the PID controller were selected. 
10 Two types of trajectories shown in Fig. 3 were applied 

as Xd(t) . One was a step trajectory having initial errors and 
the other was a sinusoidal type trajectory. 

Case 1: Step trajectory 
15 The step trajectory of equation 26 was applied to each 

axis . 



(26) 
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The control results based on the values of the PID 
controller designed according to the method of the present 
invention are given in Figs. 6A and 6B, in which response 
5 x(t) (Fig. 6A) and input torque x (Fig. 6B) are plotted 
versus time. In the upper graphs of Figs. 6A and SB, axes 1 
and 4 are represented by solid lines, axes 2 and 5 by dotted 
lines, and axes 3 and 6 by dashed lines. Likewise, the lower 
graphs of Figs. 6A and 6B have solid lines for axes 1 and 4, 

10 dotted lines for axes 2 and 5, and dashed lines for axes 3 
and 6. In all of the graphs, only solid lines are shown 
since the lines overlap. As apparent from the plots, all of 
the six axes were satisfactorily controlled when they were 
under the rule of the PID controller designed according to 

15 the present invention. 



Case 2: Sinusoidal type trajectory 

The sinusoidal trajectory of equation 27 was applied to 
each axis. 

20 

= ^rfs = ^rf6 = 240^ • sinCw^O • (l - e{- w^t)) 



The control results based on the values of the PID 
controller designed according to the method of the present 
25 invention are given in Figs. 7A and 7B, in which the response 
error e(t) (Fig. 7A) and input torque x (Fig. 7B) are plotted 
versus time. In the upper graphs of Figs. 7A and 7B, axes 1 
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and 4 are represented by solid lines, axes 2 and 5 by dotted 
lines, and axes 3 and 6 by dashed lines. Likewise, the lower 
graphs of Figs. 7A and 7B have solid lines for axes 1 and 4, 
dotted lines for axes 2 and 5, and dashed lines for axes 3 
5 and 6. In all of the graphs, only solid lines are seen since 
the lines overlap. As apparent from the plots, all of the 
six axes were satisfactorily controlled when they were under 
the rule of the PID controller designed according to the 
present invention. 

10 

As described hereinbefore, the present invention 
provides a gain selection method of PID controllers by which 
a great number of gains of PID controllers can be selected 
with a small number of parameters. As exemplified in the 

15 examples, the total number of PID controller gains to be 
selected in the six-axis robot manipulator is 18 . However, 
only two parameters (ai and suffice the selection of as 

many as 18 gains. 

Conventionally, the number of PID controller gains in 

20 an n-degree system totals 3n. However, the method of the 
present invention can reduce the parameters into one or a 
number less than n. In practice, when tuning a PID 
controller, all of 3n gains can be suitably tuned by tuning 
only one to n-1 gains even though every gain is not tuned one 

25 by one. 

Additionally, in comparison to conventional methods, 
the present invention can be applied to linear and nonlinear 
systems and even to MIMO systems. Furthermore, the present 

20 



invention enjoys the advantage of selecting and tuning gains 
of PID controllers only by the program run in the digital 
device 21 without additional equipment. 

The present invention has been described in an 
5 illustrative manner, and it is to be understood that the 
terminology used is intended to be in the nature of 
description rather than of limitation- Many modifications 
and variations of the present invention are possible in light 
of the above teachings. Therefore, it is to be understood 
10 that within the scope of the appended claims, the invention 
may be practiced otherwise than as specifically described . 
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